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I. INTRODUCTION 

It has long been surmised that the temperature of the solar 
corona and the consequent radiated power in thermal bremsstrah- 
lung is maintained by the input of mechanical energy flux of or- 
der 1 W/cm 2 from photospheric turbulent motion [1]. What is 
still a matter of contention is the physical mechanism by which 
this energy is delivered to and absorbed by the corona. It is also 
generally accepted that the solar magnetic field acts as conduit 
for the transmission of this energy flux to the corona. Of the 
three modes by which a magnetic field supports wave motion in 
a plasma only the shear Alfv^n wave can penetrate to the high 
corona and indeed along open field lines this activity is observed 
even in the solar wind. In the closed field line regions, however, 
the resolution of present observations is insufficient to detect any 
wave activity. The upper limit of the rms coronal turbulent fluid 
velocity is < 25 km/ sec [1,2]. The main difficulty is to explain the 
required dissipation of Alfv^n wave flux in the corona given the 
extremely low value of the plasma resistivity which, on the ba^ 
sis of conventional theory [3], would predict an absorption length 
orders of magnitude larger than the length of the coronal loop. 
This means that a wave-packet would bounce back and forth in 
the magnetic loop almost forever before it is absorbed. However, 
HoIIweg [3] has computed that the wave-packet is not completely 
reflected at the base of the loop in the denser chromosphere, but 
that a fraction of the energy leaks away by virtue of mode trans- 
formation, etc. If the loop is treated as a wave cavity, because of 
this leakage its quality factor Q h ~ 50 [3]. Thus, only if it can be 
demonstrated that in some manner, due to inhomogeneities, wave 
dissipation reduces Q to less than Ql can it be concluded that 
the waves deposit their energy in the corona and contribute to its 
heat content. 

Because of this inherent limitation in Alfv^n wave dissipa- 
tion, Parker has renounced [1] the hypothesis of coronal heating by 
Alfv^n waves altogether and put forward an alternative scenario. 
Since the coronal pressure is much smaller than the magnetic pres- 
sure i.e., /? = 8xp/B 2 ~ 10 -2 << 1, the equation governing the 
magnetic field is the so-called force-free equation; 

(V x B) x B = 0 . ( 1 ) 


Equation (1) has to be solved in the context of Fig. 1, where the 
feet of the field lines, presumably anchored in the dense photo- 
spheric plasma, suffer a quasi-random two-dimensional motion. 



Fig. 1 Schematic of solar coronal magnetic loop. 

This motion leads not only to the twisting of the field lines but to 
tangential discontinuities in the volume i.e., surfaces across which 
the direction of the field changes discontinuously and which are 
therefore the seat of singular current layers. In actuality, the width 
of these singular layers is determined by the resistivity, but the 
field gradients and local current densities are high, guaranteeing 
high local dissipation and ultimately to reconnected magnetic field 
lines. Parker claims that a balance is struck between winding the 
field by photospheric motion and the release of this energy by 
dissipation at the singular surfaces yielding an energy input of 
~ 1 W/cm 2 to the corona. 

Although such a process looks eminently plausible doubts 
have been expressed by van Ballegooijen [4] and Vekstein [5] on 
the basis of restricted models. In any case, analytical or numeri- 
cal demonstration of the evolution of such singular layers together 
with a quantitative determination of dissipation rates is neces- 
sary to validate this suggestion. There is, however, a very im- 
portant corollary to the evolution of singular layers viz., that the 
magnetic field can no more be considered to be regular over the 
entire volume. In the vicinity of the singular layers, neighbor- 
ing field lines will diverge exponentially over a scale length A and 
diffuse over distances greater than A with a diffusion coefficient 
Dm = ((Ax_l ) 2 /As) where As is the correlation distance along 
the field line and Ax_l is the random step in the perpendicular 
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direction. Thus the photospheric motion will give rise to irregu- 
larities in the magnetic structure of the coronal loop. 

Studies of wave dissipation in inhomogeneous but regular 
magnetic fields by Tataronis and Grossman [6], Kappraff and Tata- 
ronis [7], Heyvaerts and Priest [8], Ionson [9], Davila [10], Hollweg 
[11] and Grossman and Smith [12], reveal that the dissipation 
length is proportional to whereas in homogeneous fields it 
is proportional to Rm\ Rm , the magnetic Reynolds number, is 
~ 10 10 . If, however, the irregularities of the magnetic field dis- 
cussed above are taken into account, it has been shown by Sim- 
ilon and Sudan [13] that the dissipationlength is proportional to 
In Rm. The Q of the loop computed for such an irregular field can 
be much less than the Hollweg limit Ql so that Alfv^n dissipation 
could indeed lead to coronal heating! 

Our model for coronal heating by Alfv6n waves therefore de- 
pends on those photospheric motions with excitation frequencies 
below — (v^R* 1 to cause a slow twisting of the flux tubes 
which builds up magnetic stress and singular layers and ultimately 
to reconnection and irregular field lines. The photospheric motions 
with excitation frequencies in excess of are able to excite Alfv^n 
waves which dissipate rapidly in the presence of the stochastic field 
structure. Of course, this does not preclude direct dissipation of 
the stressed magnetic energy by joule dissipation in the singular 
layers. However, the rate of such dissipation, according to the 
Sweet-Parker model, will be proportional to Rj^ and may turn 
out not to be sufficient. 

In the next section we develop a set of rescaled MHD equa- 
tions similar, but not identical to, the Strauss equations [14,15] 
for analyzing both the slow evolution of the magnetic field and 
the fast time scale of the Alfv^n waves. The dissipation of Alfv4n 
waves is treated in Section IV. 


n. RESCALED MHD EQUATIONS FOR THE 
CORONA 

The equilibrium force balance is given by 

4*VP 0 = (V x B 0 ) x B 0 , (2) 

with P 0 = p 0 (x) - g f M dzp 0 (z), p 0 is the gas pressure, p 0 (z) is 
the stratified density and g is the acceleration due to gravity. Let 
(3 = 8xP 0 /Bq be a small quantity of order c « 1. Then, to 
lowest order in e, Eqn.(2) reverts to the force-free equation (I). 


From (1), setting B 0 = B 0 b we obtain 

V x Bo = <*oBo ; Bo • Vao = 0 , (3) 

V x 6 = c*o 6+ (6 x Vj_J9o)/Bo > (4a) 

6 • V x 6 = a 0 , (46) 

SxVxS = -i-V6 = ft R 7 1 = -(Vj_£o)/B 0 , (4c) 

V • 6 = — (6 • VBq)/Bq . (4 d) 


R c is the radius of curvature of the field line. Figure (1) shows a 
coronal loop of length L and minor radius a. We now order the 
scale of perturbations of this flux tube such that 


d/dx ± ~ 0(1) , (5a) 

d/dt ~ , (56) 


and 

(5c) 

ds L 

where 8 is the coordinate along the lines of force of Bo and x_l are 
the perpendicular coordinates. Thus, a/L is taken to be of order 
e*. Furthermore, a/R c is also of order c*. The perturbations of 
the equilibrium quantities are scaled as follows: 


p = p Q {z) + Sp(a i x jLi t) , (6a) 

p = e[p 0 (x) + 6p(s,x_L,t)] , (66) 

v = €*6v_L(*,X_L,t) + c6v||(«,X_L,t) + . . . , (6c) 

A = A 0 + e*6A||(s,xx,i) + e6Ax(s,x±,t) + ... . (6d) 

v is the fluid veloccity, A is the vector potential such that V- A = 0 
and B = Vx A. Then Ohm’s law for the plasma may be expressed 

“ a 

E = -V<p - — A =-vxB+]/ff , (7) 

where j = -(4x) _1 V 2 A is the current density and a is the plasma 
conductivity. Then, from Eqns. (6) said (7), to order , we obtain 

6vj L = (6 x V_lv?)/Bo > ( 8 ) 


and to order e taking the S component of (7), 

a 


7t SA * = 


dt 


+ (6vj_ • Vx) 


SA|1 = + , ? v i 5A ii > ( 9 ) 


with r? = (4 jrff)" 1 . Taking the 6 component of the curl of Eqn. 
(7) we get, to order e? , 


— 6P|| = (6B^ * Vj_)6t;|| + r)V 2 6B\i i 


( 10 ) 


with 

5B_l = Vj.«A|| x 6 . (11) 

To order ci the continuity equation becomes 

^Sp = -«vj. • Vp 0 , (12) 

and to order e the perpendicular momentum balance is given by: 

p 0 ±6v±=-V±6p + ti±xVo + (io + $i\\) 
x 6B± +8pg + Pov±V±6v± , 
and to order the parallel momentum balance is given by 
p 0 j t S ti|| = - V||fip + S • «jj. x «B j. + Spi ■ g + p 0 v\\ . (14) 

Here j 0 = <xoB 0 b t i/± and v\\ are the perpendicular and parallel 
kinematic viscosities, V*6vj_ is of order e, and Spt’g/pog is taken 
to be of order e$ . The pressure is obtained from a static equation 
obtained by taking the divergence of Eqn. (13) and recognizing 
that V • 6 vjl may be neglected because it is of higher order by 
e* than the remaining terms. Operating on Eqn. (13) with the 
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operator V ■ 5x we obtain to order e, 


V j. 



~Ts^ 6A " 

- Bf‘(fiB A .Vj.) V* 

+ • Vx) <*0 

+ Vx • x g 

+ I'xt'x’ Vi If, , 


(15) 


where we have eliminated 6}± through the relations Vj| • £j|| = 
“Vj. • £J_l and £j f | = -(4 jt) _ 1 V^A||. Equations (9), (10), (12), 
(14) and (15) advance, 6A\\, 8B\\, 6p t 6v\\ and 

-<P> respectively; £v x and 6B± are obtained from (8) and (11). 
This set of equations constitute a rescaled set of MHD equations 
suitable for the study of coronal magnetic fields. Equations (9), 
(12), and (15) form a self-consistent set to lowest order. 


ID. VERY SLOW MOTIONS OF CORONAL LOOP 
MAGNETIC FIELDS 


When the excitation frequency of the perturbations to the 
coronal loop is below wx (as discussed in Sec. I) Alfv^n waves will 
not be excited and the time evolution of the system is adiabatic. 
In this limit the d/dt operator is even smaller than order We 
are therefore justified in neglecting the inertial term on the LHS 
of Eqn.(15). The term involving g may also be dropped if the 
magnetic field is approximately collinear with g. Then the lowest 
order equations (15) and (9) furnish 

£ V -“> + (tT ' V.) VJ.M, - ^ V*„ ( 16 ) 

T. + (tst ■ Vi ) ” " if** + » • <«> 

where vx = Bo/i^po)^ is the Alfv&i velocity. Equations (16) 
and (17) form a closed set and describe the very slow motions 
of the loop. Nevertheless, we require the time scale r < Aj Jrj 
where Ajl is the perpendicular scale i.e., the magnetic Reynolds 
number R\f > 1. We notice that in the limit i/± — ► 0, Eqn. 
(16) is homologous to the two-dimensional incompressible Euler 
equation for the vorticity if the coordinate a is identified with 
time and the Eulerian fluid velocity with V^5A|| X S. Proceeding 
with this analogy we may identify the Lagrangian motion of the 
Eulerian fluid elements along the streamlines (magnetic lines of 
force). Thus, 

— xj. = (V_l£A|| X b)/B 0 (18) 

which leads to a Hamiltonian system of equations for the compo- 
nents x_l = (xi,Xa) 

dxi = dSA^ d\2 = dSA^ , , 

da B 0 d X 7 * da B 0 d X i ’ { } 

with 6A^/Bo as the Hamiltonian. It is well known that even if 
the velocity field in unsteady two-dimensional Navier-Stokes flow 
is continuous and laminar, the trajectories of the fluid elements 
may be stochastic, a process known as “chaotic advection” [16] or 
“Lagrangian turbulence” [17]. 

In our case 6A\\ is a function of a in addition to xi and xa 
and it is therefore highly likely that the chaotic motions of the feet 
of field lines which is equivalent to prescribing the injected current 
VI £Ajj at a = 0 will most likely lead to chaotic magnetic field 
lines. 


Equation (16) may be integrated numerically along a if the 
injected current 6 j\\ is specified at 8 = 0. The RHS of Eqn. (17) 
is therefore known and Eqn. (17) can be integrated now to furnish 
<p which specifies the velocity field. 

For Bo = Boz with B 0 constant and an injected current 
at a = 0 given by 

6j\\/jo = sin(2jrx/L) sin(2xy/L) 

(19a) 

+ € cos(2 xnx/L) cos(2tt my/L), 
we have plotted, in Fig. 2, what we define as the “irregularity” 
index 

A = \z~' ln{( | 6x ± {z) I 1 ) / ( | 5xx(0) | J > }, (196) 

which is a measure of how neighboring field lines 6xj_ apart diverge 
from each other according to d8xx/dz = B^" 1 (5B_l • V)£x_i_. In 
the limit z — ► oo, A is the Liapunov exponent. But in the case 
of finite length loops, z max = L. 


Mean Divergence Rate ^ 

w epsilon = 0.0 



z 


Fig. 2 Plot of A defined in Eqn. (19b) as a function of z for 6j\\ 
given by Eqn. (19a) for c = 0 and e — 0.3. 


IV. alfv£n waves in arbitrary magnetic 
FIELDS 

Formulation of the Eikonal Approximation 
In this case the fast component of photospheric motions excite 
Alfv^n waves so that d/dt ~ 0(e^). To lowest order we need only 
Eqns. (9), (12), and (15) for 5A||,6p, and As mentioned in 

Sect. I, the maximum bound of the perturbations in velocity in 
the coronal loop is of order 25 km/aec as compared to the typical 
Alfv4n velocity of 2000 km/aec . Thus we are justified in linearizing 
these equations to obtain 

Tt SA ii + Ts = ” v ^ a ii > (») 

JI&P = -«vx ■ Vx />0 , (21) 

Vx • *1* + SB ± • (22) 

+ -J- x g • V±6p + v ± v2 2 V±V> ■ 
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To simplify our discussion we take ao to be uniform and Sx g » 0. 
If the time dependence is taken as exp — iuit these equations may 
be combined to give: 




-«•— = 0 , 


(23) 


i.e., the distance at which the wave flux decreases by 

We first address a simple example in which the magnetic field 
Bo(x)z varies with x. Thus v A varies with x and from Eqn. (26’) 
we obtain k± = 8k\\v' A fv A = 8U} v’ A /v\ w ith v' = dv A /dx. 

Substituting this expression for k± in Eqn. (30) we get 

= (Z'’A/vu 1 v%) i oc R^f . (31) 


To solve this equation we employ the ballooning approximation 
developed in fusion physics [18,19]. In what follows we outline the 
analysis of Eqn. (23) by Similon and Sudan [13]. The perturbed 
fields are represented as 


Similarly, the dissipation layer width u; is given in order of mag- 
nitude by kj 1 at s ~ s d , i.e. 

W = (vva/ 3wv'a)$ « Rm * , (32) 


<p(x,t) = ! p ( s , x _ l ) exp[i'5(x) -«wt] , (24) 


where the amplitude <p is a slow function of s and x x while the 
phase function 5 represents the fast variation. In this represen- 
tation 6 • VS = k|} and k x = V X S and Eqn. (23) may be 
written as 


d_ 

ds 





— *1) V> = 0 , (25) 

(jJ / 


■—kx = — (VS) • kx - (kx ■ «)5 . (26) 

as 

with k = b ■ V6 . In Eqns. (25) and (26) we have set fcy — ► 0 as 
is required for long parallel wavelengths. If, on the other hand, 
one were to treat parallel wavelengths much shorter than the scale 
of variation of the equilibrium quantities with s, then d/ds is 
replaced by «fc|| in Eqn. (25) and tk± may revert back to V x if we 
wish to take account of perpendicular gradients in the equilibrium 
quantities. In this limit Eqn. (26) is modified to 


•4-kx = -(VS).kx-(k x -»c)S-(Vv A /« A )k||. (26') 
as 

The driving frequency w is constant and the boundary conditions 
are the given amplitude at s = 0 and outgoing wave condition for 
s > L. 

Dissipation of Alfv4n Waves in Complex Magnetic 
Field 

From Eqn. (25) it is straightforward to establish the wave 
energy equation 

§7 = I in I s b . (27) 

where we have dropped viscous dissipation; the wave energy flux 
F — Vx(«)e(«) and the wave energy density 

e(s) = -poSv 2 ± + (8x)~ l 6B^ 

= (8 x)- 1 [*i | | J h\ + *1 I I 2 ] (28) 

« (4x)- 1 ki | M|, | J , 

because <p « v a 6A\ From Eqn. (27) we obtain the relation 
between the wave flux at s and the wave flux at s = 0, 

F(s) = F(0)exp - f ds , (r}k 2 ± /2i/ A ) . (29) 

Jo 

We now define the dissipation length s d through 

f ds'(r}k±/v A ) = 1 , ( 30 ) 


and decreases as These results are identical with those of 

Tataronis and Grossmann [6], Kappraff and Tataronis [7], and 
Heyvaerts and Priest [8] arrived at by considerations of Alfv^n res- 
onances and matched asymptotic expansions about the resonance 
u> = k\\v A {x). The wave-packet picture described here [13] gives 
precisely the same results but is capable of greater generalization 
as we show next. 

If the magnetic field is irregular as discussed in Sect. II then 
neighboring field lines diverge exponentially, i.e., £(«o + £«) = 
£ (s 0 ) exp A 68 where £ is the distance between two neighboring 
lines. This exponentiation is limited to some correlation length 
after which it breaks down. On a scale greater than the correlation 
lengths the field lines have a diffusive behavior and £ varies as 
(2 D m 8)$ where Dm = {(Ax x ) 2 /As) is the field fine diffusion 
coefficient. The field line exponentiation stretches the wave packet 
unidirectionally and steepens the gradients. This gives rise to a 
wave packet that spreads over an area (2 Dm 8 ) constituted of thin 
filaments of width exp — As (see Fig. 3). Thus, 



Fig. 3 Wave Packet evolution in s along the stochastic magnetic 
field, (a) is the initial state; (b) and (c) show the exponentiation 
phase, during which the gradients increase as expfAs); (c) and (d) 
show the field line diffusion, during which the packet spreads on a 
length ( 2Dm »)*; (d) is the state beyond a d , when “microscopic” 
diffusivity is effective, and when energy is dissipated by resistivity 
(from Ref. [13]). 


k 2 ± (s) « k 2 exp2As. When this expression is substituted in Eqn. 
(30) we obtain 

8d = ^A -1 \n(2Xv A /rjk 2 ) a In R M > 


(33) 
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*k 0 1 is the perpendicular scale length of the perturbation at the 
photosphere a = 0. The width of the dissipation layer w is 

w & kl l (s d ) = (f?/2Atu)* oc R^} . (34) 

Significant energy dissipation does not take place until 8 ap- 
proaches 8d when the wave amplitudes crash to a low level. For 
b ~ 8<i the resistivity smooths out the large gradients, dissipates 
a large fraction of Alfv^n wave energy and redistributes the rest 
over a large area 2 Dj^Sd determined by magnetic diffusion (see 
Fig. 3). For s > the wave packet in which high k± compo- 
nents have diffused away now evolves as a result of magnetic field 
line diffusion. The wave energy flux will vary, from this point on, 
asE(s) = F(0)[l + 2D M s]- 1 . 

For oj « 2 v a /R, R m ~ 10 10 , av' A /v A = 4, R ~ 10 4 fcm 
one finds, for the laminar magnetic field, Wd/a « 1.6 x 10 -4 and 
8d/R » 780. On the other hand, for an irregular magnetic field 
with A R ~ 1/2, koa « 3x we get 8d/R = 37 and for A R ~ 2, 
»d/R « 10. The corresponding quality factors Q = Sd/ Aj| = 18 
and 5, respectively, with A|| ~ xR are easily less than Hollweg’s 
limiting Ql. Thus, an irregular magnetic field is very effective to 
dissipate Alfv6n waves because it generates the small scale lengths 
needed for dissipation. 

Dispersive Effects 

So far we have assumed that the wave packets propagate 
closely along the field lines i.e., there is no dispersion in the per- 
pendicular direction. This requires d 2 u)/dk\ « rj. A number of 
effects e.g. finite Larmor effects, finite pressure, gravity, and equi- 
librium currents could, in principle, contribute to the dispersion. 
These have been shown [13] to be insignificant in the solar coronal 
context. Because the Larmor radius of the ion is only <10 2 cm 
the FLR effect is negligible with fco 1 ~ 10 2 km. The additional 
corrections to Eqn. (23), due to finite pressure and gravity, are 
homogeneous in kj_. Furthermore, field line divergence aligns kj. 
to the direction of greatest contraction. Thus the spectrum con- 
tracts locally to one-dimension and hence d*u>/d\z±dk± tends to 
vanish leaving the wave packet nondispersive. Finally, in Eqn. 
(23) we observe that the equilibrium current does not affect the 
wave evolution if ce 0 is uniform. For arbitrary c* 0 the dispersion 
introduced by Vj_ao decreases rapidly at large wavenumbers so 
that as the wave packet develops finer structure it conforms more 
closely to the field lines. 

CONCLUSION 

We have shown that the slow motion of the feet of coronal 
arches leads to irregular magnetic fields and that Alfv4n waves 
propagating in the irregular magnetic structure are dissipated 
through filamentation of the wave packet that generates short 
scales necessary for efficient dissipation. 
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